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Abstract 

In the J, H, and Ks bands survey of the the Galactic Center region over an area 
of 2° X 5°, we have found many dark clouds, among which a distinguished chain of 
dark clouds can be identified with a quiescent CO cloud. The distances of the clouds 
is estimated to be 3.2-4.2 kpc, corresponding to the Norma arm by our new method 
to determine distance to dark clouds using the cumulative number of stars against 
J — Ks colors. Adopting these estimated distances, the size is about ~70 pc in length 
and the total mass of the cloud is 6xlO^M®. Three compact HII regions harbor 
in the cloud, indicating that star forming activities are going on at the cores of the 
quiescent CO cloud on the spiral arm. 
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1. Introduction 

The central region of our Galaxy suffers large interstellar extinction, and images taken 
toward the Galactic Center (GC) in the optical wavelengths are generally featureless star fields 
with few stars. Since the discovery of an infrared source corresponding to the GC by Becklin 
and Neugebauer (1968), the GC region has been observed extensively in the infrared, where 
interstellar extinction is one order of magnitude smaller than in the optical. Succeeding to 
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the larger-scale near-infrared (NIR) surveys of 2MASS (Skrutskie et al. 2006) and DENIS 
(Epchtein et al. 1994), which have given a general view of the Galactic structure, we have 
carried out a J, and Ks bands survey of the 2° x 5° region centered on the GC with 2 
mag deeper sensitivity and 2 times finer resolution than those of 2MASS. In the course of the 
survey, we found numerous dark clouds silhouetted against the dense star fields. 

Toward the GC, existence of many molecular clouds is well known from radio observa- 
tions (Morris and Serabyn 1996 and references therein). Their location in the line of sight is 
uncertain, but part of them, located on the near side of the GC, should appear as dark clouds 
in the near infrared images. However, the correspondence of near infrared dark clouds to radio 
molecular clouds has not been studied well. 

The distance to dark/molecular clouds is crucial for discussing their physical parameters 
associated with them. In radio observations, kinematic distances to molecular clouds are often 
derived from their radial velocity on the assumption that each cloud follows the circular Galactic 
rotation. However, since the Galactic rotation is perpendicular to our line of sight to the GC, 
we cannot determine the distance to the molecular clouds. Moreover, non-circular components 
such as expanding rings make situations complicated. 

Sofue (2006) proposed a method to determine the distance of CO clouds in the direction 
of GC from the gradient in the / — v diagram, and determined the distance of CO clouds toward 
the GC. This method is powerful for determining the distance to arm-like structures, which 
extend long in the I direction, but unsuited for small CO clouds because the velocity gradient 
is not determined well. 

In this paper, we report a chain of dark clouds found in our near infrared survey, which 
shows good positional and morphological agreements with a velocity component of ^^CO ( J=l- 
0) observed with the Nobeyama radio telescope (Oka et al. 1998). A new method to determine 
the distance to dark clouds from a diagram of J — Ks color and cumulative star number, and 
the resultant distance to the chain of dark clouds are presented. 

2. Observation 

2.1. NIR data 

Our near infrared survey toward the 2° (6) x 5°(/) region of the Galactic Center was 
carried out from 2002 to 2004. The observation was made with SIRIUS, a simultaneous imager 
in the J (A=1.26yum), H (A=1.63/im) and Ks (A=2.14/im) bands, covering an area 7'. 7 x 7'. 7 
with a pixel scale of 0".45 on the IRSF 1.4-m telescope at Sutherland, South Africa. The details 
of the instrument are described in Nagashima et al. (1999) and Nagayama et al. (2003). We 
observed totally 891 of 7'. 7 x 7'. 7 fields centered on the GC, with 10 times of 5 sec exposure. 

The reduction of near-infrared array images was performed with IRAF (Imaging 
Reduction & Analysis Facility). Each image, after subtraction of the average dark frame, 
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was divided by the normalized flat-field image. Then, the thermal emission pattern, the fringe 
pattern due to OH emission, and the slope pattern seen in the HAWAII arrays were subtracted 
from each frame with a median sky frame. The photometry for point sources was performed 
with the DAOPHOT package (Stetson 1987), and calibrated with the standard star 9172 in 
Persson et al. (1998). The lOa limiting magnitudes for the J, if, and Ks bands are 17.1, 16.6, 
and 15.6 mag, respectively, and stars brighter than 10 mag in all the three band are saturated. 
More details of our data reduction and photometries are described in Nishiyama et al. (2005). 

2.2. CO data 

i^CO (J=l-0, 115.271202 GHz) observations over the area -l.°5 < / < +3.°4 and -0.°6 
<h< +0.°6 were carried out from 1995 February to April, and from 1996 March to May with the 
45 m telescope of Nobeyama Radio Observatory by Oka et al. (1998). The spatial resolution 
of the telescope was 15" and the grid spacing was 30", resulting in an actual resolution on the 
map of 34". The details of the observation are described in Oka et al. (1998). We reprocessed 
the data and acquired 88 velocity channel maps covering the velocity range of Vlsr = —220 to 
+220 km s~^, integrated over successive 5 km s~^ widths. 

3. Comparison of NIR images with CO maps 

We examined positional and morphological correspondences between the dark clouds 
seen in our survey and the ^^CO emission. By comparing our NIR images with the velocity- 
channel maps of the ^^CO observation, we found that a filamentary high- intensity edge of ^^CO 
contours in the channel 48 (ch-48, Vlsr = +15 to +20 km s~^) shows good correlation to a 
chain of dark clouds. In Figs.l (a)-(c), we show the ^^CO intensity maps of ch-48 together 
with the adjacent channels (ch-47 and ch-49). A filamentary structure extends from (/, 6)~ 
(— 0.°4, — 0.°2) to (+0.°4, — 0.°5) only in the map of ch-48. Of particular interest is that the small 
velocity-width and velocity-gradient confined within single channel or a range less than 5 km 
s-\ 

For further comparison, we made a contour map of I(ch-48)-(I(ch-47)+I(ch-49))/2 su- 
perposed on the number density of stars brighter than 14.5 mag in the Ks band (Fig. 2), in 
which the shape of ^^CO contour conform with less star area in positional and morphological. 
Therefore, the ^^CO cloud only seen in ch-48 and the chain of dark clouds are presumably 
identical. 

4. Determination of the distance to the chain of dark clouds 

4.1. A new method to determine distances to dark clouds 

Determination of the distance to dark clouds is a key in order to discuss the physical 
properties such as its real size and mass. We will demonstrate that the cumulative number of 
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Fig. 1. i^co intensity maps for (a)ch-47 (10< Vlsb^ < 15 km s^^), (b)ch-48 (15< Vlsr < 20 km s^^), 
(c)ch-49 (20< Vlsr < 25 km-^) 

stars bluer than a certain J — K color, N{J — K), is a good indicator to estimate distances to 
dark clouds. Provided that interstellar extinction and reddening increase monotonously with 
distance in a line of sight, the observed J — K color of a star depends on the distance to the 
star. Therefore, the profile of N{J — K) against various J — is determined by extinction in 
the line of sight, and conversely, we can estimate the extinction distribution along the line of 
sight, from the N{J — K) profile. 

The number of stars to be observed can be predicted using the infrared model of the 
Galaxy (Wainscoat et al. 1992). We add an artificial screen with an extinction Ak^^dc at the 
distance D as a dark cloud, besides the general interstellar extinction Ak^ in the model. Here, 
we assume that the thickness of dark clouds along the line of sight is thin enough compared 
with the galactic scale; the thickness of screen is zero in our calculation. We then calculate 
N{J — Ks) to be observed with our actual limiting magnitudes of J = 16.0 and ^^5=15.0 mag. 
The detail of calculation is described in Appendix. 

Fig. 3 shows predicted N{J — Ks) toward = (—0.3°, —0.3°) against J — Ks- First, 
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Fig. 2. "'^^CO-intensity excess (I(ch-48) — (I(ch-47)+I(ch-49))/2) contours superposed on the Ks band star 
number density (A's<14.5; color map). Contour levels are linear with a step of 10 K km s^^. 

we discuss the case that the cloud is a complete screen and we cannot detect stars behind the 
cloud. In Fig. 3 (a), the 10 lines {D = 1.0 to 6.0 kpc with a step of 0.5 kpc) are shown, together 
with the line without a screen. The top thick line represents the case of no-screen with only 
interstellar extinction, while the 10 lines are for a complete screen with Aks,dc = oo in addition 
to interstellar extinction. In these cases, N{J — Ks) is contributed by only foreground stars. 

The N{J — Ks) profiles with the complete screen increase as same as the no-screen 
N{J — Ks) profile at small J — Ks because of the foreground stars, but are detached from the 
no-screen line, beginning with the nearest D case. Because the detached points are different 
for different D, we can determine D by comparing predicted and observed detached points. 
After detached from the no-screen line, the increase of N{J — Ks) stops at Nfg, the number of 
foreground stars detectable with our sensitivity. Njg is uniquely determined for D in the cases 
of completely obscuring clouds. 

Next, we discuss the case that the cloud does not obscure the stars behind completely. 
In this case, certain number of stars behind the cloud are detected with extra extinction and 
reddening due to the cloud. Fig. 3 (b) shows the case of -D = 2 kpc with Aks,dc = 0-1 — 1-0 mag 
(0.1 mag step), 1.2, 1.5 and oo mag. As we described in Fig 3 (a), N{J — Ks) with Aks,dc = oo 
stays at Nfg because background stars are undetectable. In the case of Aks,dc=^-0 mag, all 
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J-Ks [mag] 
(b) the case for non-complete obscuration 

Fig. 3. The predicted N{J — Ks) versus J — Kg color by changing screen parameters; (a) distance D, 
from 1.0 to 6.0 kpc (0.5 kpc step), and (b) extinction Ak,dc from 0.1 to 1.0 mag (0.1 mag step), 1.2, 1.5 
and oo. The thick (upper most) hnes in each panel indicates the N[J — Ks) without a screen (interstellar 
extinction only). 

background stars are about 2 mag {Ak^ = 0A94:E{J — Ks), Nishiyama et al. 2006a) reddened 
in J — Ks by the screen. Therefore, the increase of the N{J — Ks) with Aks,dc=^-0 stops 
at Nfg at J — Ks=l, but N{J — Ks) increases again from J — Ks=S due to contribution of 
background stars. Because the different Aks,dc gives different width of J — Ks stopping at 
Nfg, we can determine A^g from it. 

As we described using Fig. 3 (a) and (b), the N{J — Ks) profiles reflects D and Aks,dc, 
separately. Therefore, we can determine the D and Axg^DC by comparing the observed N{J — 
Ks) with the template parameterized by D and Aks,dc- Fig. 4 shows simplified templates to 
determine D and Aks,dc for {h^) = (—0.3°, —0.3°); a total of 13 lines for the cases of no dark 
screen (interstellar extinction only) and combination of Ak^ jjc (0.5 and 1.0 mag) and D (1, 2, 
3, 4, 5 and 6 kpc) are shown. Because the lines with various combination of Axg.nc and D are 
well separated from each other, we can assign an unique pair of Axg^DC and D for the observed 
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Fig. 4. A simplified temperate to determine D and Aks,dc using predicted N{J — Ks) are shown. The 
hnes with A^^ =0.5 and 1.0 mag located at £'=l(red), 2(orange), 3(green), 4(light green), 5(blue), and 
6(light blue) kpc are drawn. Upper and lower tracks of each color indicate Aks,dc~^-^ aud 1.0 mag, 
respectively. 

N{J — Ks) profile of dark clouds. 

4-2. Comparison of our method with the Wolf diagram 

In the Wolf diagram, a traditional method to estimate the distance to dark clouds, 
cumulative number of stars observed brighter than magnitude m, N{m), is adopted instead of 
N{J — K). The distance to the dark cloud is estimated from the apparent magnitude where 
the N{m) profile for the dark cloud is detached from that for the comparison field. In reality, 
however, the detached point is very ambiguous due to the large dispersion in the absolute 
magnitudes of the stars as pointed out by Trumpler and Weaver (1953). 

By contrast, the intrinsic color of stars in the near infrared tends to be similar among 
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almost all the spectral types; the J — K of dwarfs and giants lie between 0.0 and 1.0 (Wainscoat 
et al., 1992). The color excess E{J — K) due to interstellar reddening is also expected up to ~ 4 
mag at the distance of the GC. Thus, the large change oi J — K compared with the ambiguity 
of intrinsic J — K oi stars enables us to determine the distance to dark clouds. 

4-3. Distance to the chain of dark clouds 

Ten zones were selected around and in the chain of dark clouds : two reference fields (RFl 
and RF2) and eight dark clouds (DCl to DCS) (Fig. 5). In order to confirm the reliability of 
the Galaxy model toward the GG region, we compared the observed and calculated N{J — Ks) 
toward the reference fields, RFl and RF2. The N{J — Ks) profiles for two reference fields are 
both in fairly good agreements (Figs. 6 (a) and (b)). We therefore regarded the Galaxy model 
as reliable enough and adopt it for determination of the distances. 

We calculated N{J — Ks) for DGl to DG8. We tried to search for the best fit variables, 
^Ks,DC and with steps of 0.2 mag and 0.1 kpc, respectively. The model lines with best fit 
Aks,dc and D for DGl-8 are shown together with the observed N{J — Ks) histograms (Figs. 
6 (c)-(j) ). The best fit lines show good agreement compared to the RFl and RF2, where the 
model lines show small discrepancies. A possible reason that the fittings for RFl and RF2 are 
worse than those for the dark cloud regions is the uncertainty of model for the inner-galaxy, 
which becomes insignificant when dark clouds are present. 

The results for DG1-DG8, summarized in Table 1, show convergence to a specific dis- 
tance of 3.2-4.2 kpc. Presumably, these dark clouds lie on a single chain of molecular clouds at 
a distance of ~4 kpc with the extinctions range of Aks,dc ~ 0.4 to 1.0 mag. The interstellar 
extinction in front of the clouds, Aksi is estimated to be 0.4 to 0.5 mag from the model calcu- 
lation. Most of the dark clouds previously studied are nearby, and the foreground extinction 
is negligibly small. In contrast, the dark clouds studied in this paper are significantly farther, 
and therefore the cumulative interstellar extinction is comparable to Akq^dc- 

The GO cloud discussed in this paper was also identified with Vlsr =21 km s~^ and 
dv/dl = 5.1 km s~^/deg (Sofue 2006). Its distance was derived to be 4.7 kpc for Rq =8.0 kpc, 
which is now 4.4 kpc for Rq =7.5 kpc, and is slightly larger value than ours. 

The distance 3.2 — 4.2 kpc is halfway to the GG at the distance of 7.5 kpc from us. 
According to Vallee (2008), a summary of previous studies of galactic structure, three arms, 
Sagittarius-Garina, Scutum-Grux and Norma arms intervene between the Galactic Genter and 
us. In Fig. 1 of Vallee (2008), these three arms are placed to cross with the line of sight 
to the GG at approximately 0.8, 2.8 and 4.0 kpc from us, respectively. It is reasonable that 
the dark/molecular cloud is located around Norma arm, and rotates around the GG with a 
non-circular velocity components Vlsr ~ +20 km s~^ toward the GG. 

The Four-kpc molecular ring (Glemens, Sanders and Scoville 1988, Nakanishi and Sofue 
2006) is also known as a component, crossing the line of sight to the GG. However, Vallee (2008) 
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Fig. 5. Positions of RFl-2 and DCl-8 in the galactic coordinate, together with CO contour. 

claimed that the Four-kpc molecular ring is just assemblage of starting segments of four spiral 
arms, and we therefore leave it out of consideration. 

5. Discussion 



We estimated a projected length of the chain of dark clouds to be ~70 pc and the typical 
thickness of the filaments is as thin as 0.1 degree, or about 7 pc, located at ~30 pc off from the 
Galactic plane. The gaseous filaments are themselves divided into finer structures, and several 
high-intensity knots are seen on the ridge of the ^^CO contours. By integrating the ^^CO excess 
in the Fig. 2, the masses of each ridges are estimated to be an order of ~10^M ® , and the 
total mass for the entire CO cloud results in 6xlO^M®. If we assume the same thickness in 
the line of sight (7 pc) as the apparent thickness and an extent of 70 pc, this total mass yields 
molecular gas density on the order of 18 M®pc~^, or 4x10^ H2 cm~^. This is a typical density 
for normal molecular clouds. 

In the three high-intensity knots, we found four NIR sources, SFl-A, SFl-B, SF2, and 
SF3, with extents of ~ 10" to ~ 40" or 0.2-0.8 pc at 4 kpc. Pseudo color images and position of 
them are shown in Fig. 7 and Table 2. Mid infrared sources (MSX-PSC, Egan & Price 1996) 
are associated with all of them. OH, H2O, and/or CH3OH masers (SIMBAD, or Caswell & 
Haynes 1983, Caswell et al. 1983 and Menten 1991) are also associated with SFl-A, SFl-B, 
and SF3 with Vls'r=+15 to +23 km s~^, quite consistent to Vlsr of ^^CO. Both their sizes and 
Ks magnitudes in 2MASS Point Source Catalog (9.15, 7.67, and 7.80 for SFl-A, SF2, and SF3) 
are consistent with those of typical compact Hll regions at a distance ~4 kpc. We infer that the 
dark/molecular cloud discovered in this survey are molecular clouds in the Norma spiral arm 
at 4 kpc from us. 
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Fig. 6. Cumulative star-number histogram against J — Ks toward two reference fields(RFl and RF2) and 
eight DCs (DCl-8) arc drawn. The dashed line indicates the predicted number from the Galaxy model. 



6. Summary 



We have found a number of dark clouds in our near infrared survey in the direction of 
the Galactic center region over an area of 2° x 5°. From the agreement between the infrared 
and CO morphologies, we have identified a chain of the dark clouds located around the Norma 
spiral arm with a single CO molecular cloud lying at the distance of ~4.0 kpc. The chain of 
dark cloud harbors compact HIT regions at the cores of the cloud in the CO spiral arm. 
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Appendix 1. The Galaxy model: Stellar populations with interstellar extinction 

The numbers of stars detectable in the J, and Ks bands have been estimated on the 
basis of the Galaxy model of Wainscoat et al. (1992). It comprises geometrically and physically 
realistic representation of the Galactic disk, bulge, halo, spiral arms, and a molecular ring. They 
represent each of the distinct Galactic components by space densities of 87 spectral-type stars 
with scale heights and absolute magnitudes at B,V,J,H,K, 12 and 25 fim. 

In our calculation, three components, an exponential disk, a bulge, and a molecular 
ring, were taken into account. We adopted 7.5 kpc as Rq, the distance between the GC and 
us (Eisenhauer et al. 2005; Nishiyama et al. 2006b) instead of 8.5 kpc used in Wainscoat et 
al.(1992), which leads to a reduction factor 0.88 for the associated parameters. We also took 
only 29 spectral types of main sequences and giants from Table 2 of Wainscoat et al.(1992) into 
calculation. 

The stellar flux is subject to general interstellar extinction. The amount of general 
interstellar extinction between us and distance Dobj, ^Dobji calculated as 



Here, r and z are the distance from the Galactic Center on the Galactic Plane and the distance 




(Al) 
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Table 1. Positions of DCl-8 and RFl-2 



ID 


I (deg) 


b (deg) D (kpc) 


Aks,dc (mag) 


DCl 


+0.250 


-0.450 


4.2 


1.0 


DC2 


+0.100 


-0.465 


4.2 


0.6 


DCS 


-0.025 


-0.445 


3.2 


0.8 


DC4 (SF3) 


-0.200 


-0.445 


3.8 


0.6 


DC5 (SF2) 


-0.290 


-0.400 


4.2 


0.4 


DC6 (SFl) 


-0.395 


-0.225 


3.6 


0.8 


DC7 


-0.095 


-0.305 


4.0 


0.8 


DCS 


+0.185 


-0.315 


4.2 


0.6 


RFl 


-0.230 


-0.270 






RF2 


+0.250 


-0.650 






Table 2. Positions of SFl- A, 


SFl-B, SF2, SF3 




ID 


I (deg) 


b (deg) 






SFl-A 


359.615 


-0.244 






SFl-B 


359.654 


-0.241 






SF2 


359.716 


-0.374 






SF3 


359.370 


-0.456 





from the Galactic Plane, and hr and hz are scale length of r and z, adopting 3080 pc and 100 
pc, respectively, jgc is an amount of extinction per unit length at the Galactic Center, and it 
was obtained as follows. 

First, we obtained the averaged extinction of red clump stars in the RFl and RF2 
(y4^g=1.7 and 1.2, respectively), and assumed them as the total amount of general interstellar 
extinction to the Galactic Center [D obj =7. 5kpc). The detail of method to obtain extinction is 
described in Nishiyama et al.(2005). Next, we substituted Ao^^j = 1-7 for RFl and 1.2 for RF2 
into the above equation, and derived 7gc=0.76 for RFl and 0.71 for RF2. We adopted 0.73 
mag/kpc, the average of the above two values, as 7gc in this work. 

We calculated the apparent magnitude of 29 spectral types of stars every 10 pc bin 
along the given line of sight, and integrated the star number detectable by the given limiting 
magnitude and color (J < 16 and Ks < 15, 85% detection rate of artificial stars in the most 
crowded region of our observation) over the distance of 20 kpc. We consider their sum to be 
the total number of stars to be detected. 
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(a) SF i-A and i-B 



(c) SF 3 



Fig. 7. Composite color (J: blue, H: green, and Ks: red) images of SFl, 2, and 3. North is up, east is 
left. Scale is the same among them. 
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